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While the phase of a coherent light field can 
be precisely known, the phase of the individual 
photons that create this field, considered individ- 
ually, cannot Phase changes within single- 
photon wave packets, however, have observable 
effects. In fact, actively controlling the phase of 
individual photons has been identified as a pow- 
erful resource for quantum communication pro- 
tocols [2|, l3|] . Here we demonstrate the arbitrary 
phase control of a single photon. The phase mod- 
ulation is applied without affecting the photon's 
amplitude profile and is verified via a two-photon 
quantum interference measurement 0? S|? which 
can result in the fermionic spatial behaviour of 
photon pairs. Combined with previously demon- 
strated control of a single photon's amplitude 
P, jzL [H, @, Eoj], frequency fllj ]. and polarisation 
[12|, the fully deterministic phase shaping pre- 
sented here allows for the complete control of 
single-photon wave packets. 

Consider two identical photons mode-matched at the 
two input ports (A and B) of a 50/50 non-polarising 
beam splitter (NPBS), represented by the initial state 
= ( see Fig. [T]). Due to the indistinguishabil- 

ity of the photons, the detection of one photon in out- 
put port C or D at time to projects the initial prod- 
uct state into the "which path" superposition state 
|*±(t )> = (|U,0b> ± \0aAb))/V2 of the remaining 
photon. As first demonstrated by Hong, Ou and Man- 
del the bosonic nature of photons always results in 
the detection of the second photon in the same output 
port as the first. However, we can alter this coalescence 
behaviour by introducing an arbitrary differential phase 
between the two components of |\P±)- This results in 
a phase-dependent wave function of the remaining single 
photon 

|*±(*o + r)) = (\1 A , Ob) ± e iA ^\0 A , \ B ))/y/2. (1) 

From (pQ) , it can be readily shown that the probability to 
detect the two photons in opposite output ports depends 
on the magnitude of the applied phase shift between the 
detection times 
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Pcoinc = sin 2 (A^/2). 



(2) 



Hence, with well tailored phase shifts, we gain control 
over the statistics of coincident photon detections. 

In our experiment, single photons are generated from 
a coupled single-atom cavity QED system (see Fig. [T] and 
the Methods section). This apparatus allows us to pro- 
duce photons at the output mode of the cavity with a 
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FIG. 1: Experimental apparatus. A single 87 Rb atom 
is coupled to the optical mode of a high-finesse cavity, which 
when subjected to a sequence of laser pulses, emits a stream 
of single photons. A polarising beam splitter (PBS) directs 
the photons randomly into one of two optical paths. In the 
first path, an electro-optic modulator (EOM) is used to arbi- 
trarily shape the phase of the single photon wave packets. A 
delay fibre in the second optical path allows two subsequently 
generated photons to synchronously impinge on a 50/50 non- 
polarising beam splitter (NPBS). The photons are detected 
by single photon counting modules (SPCMs) which are con- 
nected to a time to digital converter (TDC). 



well defined spatio-temporal mode profile @, 0, H, 0] • In 
particular, we design the photons with a Gaussian tem- 
poral profile — a shape proven to be optimal for two 
photons interacting on a beam splitter [141 ] . Our pho- 
ton wave packets have a full width at half maximum of 
150 ns, much longer than our single photon detection 
time resolution of 2 ns. The single-photon character of 
our photon stream is verified by measuring the second- 
order correlation function g^ 2 \r). The resulting suppres- 
sion of coincidence events of 95 ± 1% is consistent with 
independently measured influences of multi-atom events 
within the cavity and noise on the single photon counting 
modules (SPCMs). 

The single photons are generated with a random po- 
larisation such that a polarisation beam splitter (PBS) 
equally distributes them into two optical paths (Fig. [T]). 
In one path, a fibre based electro-optic phase modula- 
tor (EOM) is used to arbitrarily shape the phase of the 
photonic wave packets. The other optical path consists 
of a 270 m long delay fibre. These two optical paths 
are finally mode-matched on a NPBS, after which the 
photons are detected by SPCMs connected to a time-to- 
digital converter. The photons are generated at a rate 
of 740 kHz, chosen to match the path length difference 
between the EOM fibre and the delay fibre. With this 
configuration, two independent, subsequently generated 
Dhotons can arrive at the beam splitter simultaneously 



pnotons cai 

[H ElIEDI. 



In a first experiment, we demonstrate the ability to 
convert the usual bosonic coalescence for photon pairs 
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FIG. 2: Photon pair coincidences, a, One of the two 

photons of each pair is subjected to a 7r-phase shift at the 
centre of the wave packet envelope via a voltage step applied 
to an EOM. b,c, Connected dots refer to measured coinci- 
dences for non- interfering photons (see the Methods section) 
and histograms refer to measured coincidences of interfering 
photons, b, Interfering photons detected in the same tempo- 
ral half of the photon wave packets (both within I or both 
within II) should never be detected in opposite output ports. 
The total measured coincidence detection rate is only 16 ±2% 
of that for non- interfering photons, c, Interfering photons de- 
tected in opposite halves of the wave packet (one within I and 
the other within II) are ideally twice as likely to be detected 
in opposite output ports as compared to non- interfering pho- 
tons (i.e. always). The measured coincidence rate is 183 ±6% 
that for non-interfering photons. All stated uncertainties are 
statistical. Due to the finite rise time of the phase shift, pho- 
tons detected within ±5 ns of the step are omitted from the 
above data. Photons detected further than 80 ns away from 
the step are also omitted due to the low signal to noise ratio. 



into anti-coalescence by introducing a relative phase shift 
of A^ = 7r between the photon detections. This is ac- 
complished by subjecting the photon traveling through 
the EOM to a sudden tt phase shift exactly at the centre 
of the wave packet envelope (Fig. [2]). This results in the 
total wave functions of the two photons to be orthogonal 
to each other. Hence, when integrating over the entire 
wave packets, it is equally likely to detect the two pho- 
tons in the same or opposite detectors. However, closer 
examination reveals that the photon counting statistics 
depend on when the photons are detected with respect to 
the applied phase shift, as only a phase shift between pho- 
ton detections is relevant. For instance, if both photons 
are detected within the same temporal half of the wave 
packet envelopes (as defined in Fig.E^i), then there is no 
relative applied phase shift (A^ = 0), and the photons 
coalesce. On the other hand, if the photons are detected 
in opposite temporal halves (A^ = 7r), then according to 
equation ([2]), the photons are always detected in oppo- 
site output ports. These coincidences are therefore twice 
as likely as for non-interfering photons, equivalent to a 



200% 



150% 



100% 




50% 



EOM voltage change (V) 

FIG. 3: Coincidence rate vs. applied phase shift. Red 

dots show the coincidence detection probability with respect 
to non- interfering photons for photon pairs detected in op- 
posite temporal halves of the photon wave packet. The solid 
fitted line follows the sinusoidal behaviour of equation J2} with 
a visibility of 67 =b 3%. The coincidence rate for photon pairs 
detected in the same half of the photon wave packet remains 
at a low value of « 14% (individual points: blue squares; av- 
erage value: solid blue line). An offset of 5% due solely to 
multi-atom effects and noise on the SPCMs is indicated by 
the filled area. The phase shift of the EOM is proportional 
to the applied voltage, verfied via an independent measure- 
ment with continuous laser light, with a resulting tt- voltage 
of 3.60 =b 0.05V. The data from Figure ^jp-c is incorporated 
into this Figure. 



fermionic spatial behaviour. 

The results of our experiment are shown in Figure [2Jd- 
c. When both photons are detected in the same temporal 
half of their wave packets, they exhibit a strong coales- 
cence resulting in a coincidence rate only 16 ±2% of that 
for non-interfering photons (see the Methods section). 
On the contrary, photons detected in opposite tempo- 
ral halves show anti-coalescence, where the detection of 
photon pairs in different output ports is 183 ±6% as likely 
as for non-interfering photons (Fig. [2k). The deviation 
from the respective ideal values of 0% and 200% can be 
explained by our non-perfect single photon source as well 
as by mode mismatch on the beam splitter. In the lat- 
ter case, the spatial and polarisation modes overlap with 
greater than 99% fidelity. However, shot-to-shot varia- 
tions in the amplitude and frequency profile of the emit- 
ted photons, due largely to addressing different Zeeman 
sub- levels in the photon generation process, can lead to 
a temporal or frequency mismatch on the beam splitter 



As mentioned previously, when integrated over the en- 
tire wave packets (sum of the data from Figs. [2)3 and c), 
the interference effects vanish and the total detected co- 
incidences are 98 ± 3% compared to non-interfering pho- 
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tons. These results are analogous to the case of orthogo- 
nally polarised photons — while a vertically and a hori- 
zontally polarised photon are non-interfering on a NPBS, 
perfect coalescence (anti-coalescence) will still arise when 
the photons are detected with parallel (orthogonal) po- 
larisations in diagonal bases [161 ]. 

Next, we extend our experiment to investigate the ef- 
fect of discrete phase shifts with arbitrary magnitude 
A«£. As seen in Figure 03 photons detected in the same 
temporal half of the wave packet envelopes always coa- 
lesce, independent of the applied phase shift. However, 
photons detected in opposite temporal halves exhibit a 
sin 2 (A0/2) dependence of the photon pair coincidence 
detection rate, as given by equation ([2]), with a visibility 
of 67 ± 3%. As expected, a maximum number of coinci- 
dences is observed for phase shifts of odd multiples of 7r, 
whereas the photons coalesce for even multiples of tt. The 
fact that our visibility remains consistent even for high 
applied phase shifts testifies to the stability of our phase 
modulation, determined mainly by the voltage generator. 
In the case of = 0, 27r, the coincidence rates for pho- 
tons detected in opposite temporal halves of the photon 
wave packets should be the same as for those detected 
in the same temporal half. However, we observe a small 
discrepancy due to accumulated phase noise from the fre- 
quency and amplitude uncertainty in our single photon 
stream. This leads to an increased coincident detection 
rate for hi gher detection time differences, as discussed in 
references (ill [l3| . 

In addition to discrete phase shifts, we also investigate 
the general case of continuous phase shaping within a 
photon wave packet. As an example, applying a linear 
phase ramp to the photon traveling through the EOM 
results in a frequency shift, where the magnitude of the 
shift is controlled via the slope of the ramp. For two 
photons with a frequency difference Ah> detected at a time 
difference r, the relative phase shift is A^(r) = 2ttAi/t. 
Hence, the probability of photon coalescence oscillates 
sinusoidally at the difference frequency as a function of 
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the detection time difference 

We demonstrate this phase shaping with two different 
phase ramp experiments. In the first experiment, we ap- 
ply a single phase ramp from to 4.47T over the central 
200 ns of the photon in the EOM fibre path, resulting 
in a frequency shift of 11 MHz (Fig. [4^). In the second 
experiment, making use of the periodicity of equation 
([2]), a frequency shift of 25 MHz is accomplished via a 
saw-tooth function where each step of 40 ns introduces a 
phase increase of exactly 2tt (Fig. HJd). While the finite 
fall time of the saw-tooth function may potentially lead 
to phase errors, this technique allows for arbitrarily high 
frequency shifts by requiring an absolute phase range of 
only 27r on the EOM. In both experiments, the measured 
signal for interfering photons shows the expected oscilla- 
tions with high visibility (Fig. [4]>d). 

The presented phase modulation scheme not only re- 
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FIG. 4: Linear phase ramp. A linear phase ramp produces 
a frequency shift of the single photons, verified via a quantum 
beat with a second, unmodulated photon, a, A frequency shift 
of 11 MHz is applied via a linear phase ramp of 4.47T over the 
central 200 ns of the generated photon wave packets (gray 
shaded curve), b, A frequency shift of 25 MHz is applied by a 
series of 40 ns long 2tt linear phase ramps. c,d, Open red dots 
show the measured coincidences for interfering photons sepa- 
rated in frequency by 11 MHz (c) and 25 MHz (d). Closed 
black dots show the reference data for non-interfering pho- 
tons. The black solid lines show the Gaussian shapes for non- 
interfering photons, calculated from the width of the photon 
wave packets and using only the amplitude as a fit parame- 
ter. The red solid lines are the predicted curves for interfering, 
frequency-shifted photons, calculated directly from twice the 
amplitude of the result for non- interfering photons modulated 
by equation J2} with the appropriate frequency. All photons 
detected within the 300 ns windows shown in (a) and (b) are 
included in the data. 



veals fundamental aspects of two-photon interference, 
but also has various potential applications. For instance, 
shaping the phase of the reference photon would al- 
low one to fully determine the phase properties of an 
unknown photon using this technique. Similarly, this 
method to measure the phase change of single photons 
can be utilised for quantum key distribution protocols 
[3, 0] , in particular for time-bin entanglement [l8[ , with 
the advantage of not requiring interferometric stability. 
Finally, as any phase shape can be created or even com- 
pensated, this technique may be utilised fo r q uantum 
computation protocols with linear optics 19|, |2(| . 



METHODS 

Cavity QED based single photon source. Indi- 
vidual 87 Rb atoms are probabilistically transferred from 
a magneto-optical trap into a high-finesse optical cavity 
(F « 6 x 10 4 ) via a horizontally aligned dipole trap guid- 
ing beam at 1064 nm oriented orthogonal to the cavity 
axis . The transit time of an atom through the cavity 
mode is « 140 /is, and the probability to have one atom 
in the cavity at a given time is about 3%. 

To generate single photons, the intracavity atom is op- 



tically pumped to the 2 Si/ 2 F = 2 hyperfine ground 
state via a 800 ns long laser pulse resonant with the 



S1/2 F = 1 <-> P3/2 F' = 2 transition. Next, a laser 
pulse resonant with the F = 2 <-> F' = 1 transition to- 
gether with the cavity resonant with the F = 1 <-> F' = 1 
transition drives a vacuum-stimulated Raman adiabatic 
passage 21, 22| resulting in a single photon in the cavity 
mode and the atom in the F = 1 state. The cavity is 
optically asymmetric (mirror transmissions of 2 ppm and 
101 ppm) such that the photon is preferentially emitted 
through the higher transmission mirror into the spatial 
output mode defined by the cavity. The amplitude pro- 
file of the emitted photon is controlled by the temporal 
shape of the F = 2 <-> F' = 1 excitation laser pulse. 
With an atom present in the cavity, the probability of 
generating a single photon in a given trial is « 7%, and 
the overall transmission and detection probability of this 
photon is also ~ 7%. The insertion loss of our particular 
EOM is about 5 dB. 

Non-interfering photons. In this manuscript, the 
term non-interfering photons refers to photon pairs ar- 
riving at the NPBS with no temporal overlap or with 
orthogonal polarisations. Such photons exit the NPBS 
independently and are therefore equally likely to be de- 
tected in the same or opposite output ports. For our 
reference curves, we use the average number of coinci- 
dences in the neighbouring peaks (2 < \n\ < 6) in the 
second-order correlation function. With equal intensities 
at both input ports of the NPBS (measured to be accu- 
rate within ±5%), these coincidences occur twice as often 
as coincidences within the n = peak and therefore are 
divided by 2 for comparison [23[. For \n\ = 1, the num- 
ber of correlations is 75% that of the outer peaks due to 
the probabilistic coupling of photons into either the long 
fibre or the EOM fibre by the PBS. 
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